On the absence of the glass transition in two dimensional hard disks 



O 

o 

(N 

m 
o 



c3 



^3 

o 
o 



> 

OS 

o 
o 

o 

-i— > 

c3 



i 

c 

o 
o 



X 



Marco Tarzia a ' b 

Dipartimento di Scienze Fisiche and INFN sezione di Napoli, Universitd degli Studi di Napoli "Federico II", 

Complesso Universitario di Monte Sant'Angelo, via Cinthia, 80126 Napoli, Italy and 
b School of Physics and Astronomy, University of Manchester, Manchester, M13 9PL, United Kingdom 

(Dated: February 6, 2008) 

In this paper we study the glass transition in a model of identical hard spheres, focusing on the 
two dimensional case. In the mean-field limit the model exhibits an ideal glass transition of the same 
nature of that found in discontinuous spin glasses. Nevertheless, a systematic expansion around the 
mean-field solution seems to indicate that the glass transitions is smeared out in two dimensions, 
in agreement with some recent results. Our investigation could be generalized to higher spatial 
dimensions, providing a way to determine the lower critical dimensionality of the mean-field ideal 
glass picture. 
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Under fast enough cooling or densification diverse ma- 
terials, such as molecular and polymeric liquids, col- 
loidal suspensions, granular assemblies, molten mixtures 
of metallic atoms, may form glasses i.e., amorphous 
states that may be characterized mechanically as a solid, 
but lack of the long range crystalline order. Despite all 
the work devoted to the subject, the underlining mech- 
anisms responsible for the vitrification processes are not 
well understood, as the transition to the glassy state is 
still deemed specifically to be one of the most obscure 
enigmas in condensed matter physics. Many valuable 
theories attempt to describe these remarkable phenom- 
ena, but none of them is as yet regarded as compelling. 

A system of identical hard spheres confined in a fixed 
volume [1, d, 0, H, [1, 0, IE is the simplest system ex- 
hibiting a dramatic slowing down of the dynamics in the 
high volume fraction region, referred by many authors as 
a glass transition [H, |3,H, @| • In a recent paper [§[ , us- 
ing a replica approach it has been very nicely shown that 
hard spheres in large space dimensions undergo an ideal 
glass transition at a volume fraction 4>k ■ The same result 
has been obtained in three dimensions employing diverse 
kinds of mean-field-like approximations (such as the hy- 
pernetted chain approximation Q and the small cage ex- 
pansion 0,i) and using the Carnhan-Starling equation 
of state. However, for finite dimensional systems (with 
short range interactions) the mean-field picture should be 
modified by (non-perturbative) activated events [l(| , and 
one might wonder to what extent the mean-field scenario 
is still valid in that case. 

Here we study a system of identical hard disks on a 
two dimensional square lattice. We first analyze the 
mean-field solution, where an ideal glass transition of 
the same nature of that found in mean-field model for 
glasses occurs. We then consider a systematic ex- 
pansion around the mean-field limit, which allows to take 
into account short range correlations as corrections to the 
mean-field approximation. This is accomplished by con- 
sidering bigger two dimensional L x L square cells of size 
L = 2,3,..., within which the model is solved exactly 
(see also Ref. [H, EH)- We observe that the glass tran- 



sition occurs at higher densities as the size of the cell 
is increased, and seems to be smeared out in the limit 
L — > ao. This analysis hints that there is no glass transi- 
tion in two dimensional hard disks, in agreement with the 
recent results of Refs. [aly, |7j. Interestingly enough, the 
study presented here could be easily generalized to higher 
spatial dimensions, providing a direct way to estimate the 
lower critical dimensionality of the mean-field ideal glass 
picture. The latter investigation could prune down the 
number of candidate theories for the glass transition. 
The Hamiltonian of the model reads: 



H 



^ J riitij, 
[ij] 



(1) 



where the lattice variables rii = 0, 1 whether the cell i is 
occupied by a particle or not, and the sum is restricted 
over the couples of sites such that their distance is 
equal or less than two lattice spacings: dij < 2. The 
model can be regarded as a system of hard disks of di- 
ameter y/b—e lattice spacings. The limit J — > oo is taken, 
insuring the hard core exclusion. 

The model can be solved in mean-field on the random 
regular graph [l4j, i.e., a random lattice where every ver- 
tex has k + 1 neighbors, but which is otherwise random 
(to mimic the two dimensional case, we choose k = 3). 
Locally the graph has a tree-like structure with a finite 
branching ratio, but has loops of typical size IniV. The 
presence of loops is crucial to insure the geometric frus- 
tration, but the local tree structure allows for an ana- 
lytical solution of the model, since we can write down 
iterative equations on the local probability measure. To 
this aim, let us consider a branch of the tree ending on 
the site i, and denote by ij, j G {1, . . . , A;}, the roots of 

the branches connected to the site i. We call the 
partition function of this branch restricted to the config- 
urations where the site i is occupied by a particle. Analo- 
gously, we define the partition function of the branch 
restricted to configurations where the site i is empty, and 

Z the partition function of the branch restricted to con- 
figurations where the site i is empty with all neighbors 
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FIG. 1: Mean field solution of the model on the random graph. 
Density, p, and entropy per site, S, as a function of the chem- 
ical potential, fj,, in the liquid (black continuous curve), crys- 
talline (blue dashed curve), and glassy (red circles) phases. 



ij empty. Using the Grand Canonical ensemble, the fol- 
lowing recursion relations are derived: 

k 
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where p is the chemical potential. It is convenient to 
introduce on any site i the local cavity fields (3hi = 

]n{Z^/Z^) and 0Vi = ln(Zp ) /Z^), in term of which 
the iteration relations, Eqs. @, read: 
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3 = 1 



Ph< 
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Starting from Eqs. ([3]), we find at low density a liquid 
phase, characterized by a homogeneous (replica symmet- 
ric) solution, hi = h and u,; = v. Given this solution, 
using Eqs. ^ and (0, the thermodynamic quantities 
can be derived, and, in particular, the density p = (ri^) 
and the entropy per lattice site, S = —f3F — [3pp, are 
obtained (continuous curve in Fig. QJ. As the chemical 
potential (the density) is increased, a first order phase 
transition from the liquid phase to a crystalline phase 
occurs at a melting point p m ~ 2.46 (p m ~ 0.1364). 
The crystalline state is characterized by a periodic struc- 
ture, which breaks down the translational invariance, and 
can be obtained introducing different sub-lattices. More 
precisely, we introduce three sub-lattices, a, 6, and c, 
on which the local cavity fields are site independent. 
The sub-lattices must be organized in such a way to re- 
produce the crystalline order: each vertex of the sub- 
lattice a is connected with 3 sites of the sub-lattice 6 
[i.e., a — > (6,6,6)]. Analogously, we have 6 — > (c, c, c) 
and c — > (c, c, a). Also the free energy shifts, Eqs. ^ 
and (O, must be computed carefully, taking into ac- 
count the structure of the three sub-lattices. In the 
crystalline phase, as the chemical potential is further in- 
creased above p m , the density rapidly approaches the 
maximum density, p mQX = 0.2, and the entropy per site 
approaches zero (dashed curves in Fig. [I} . 

The crystallization transition can be avoided and, in 
this case, the system enters a supercooled state, still de- 
scribed by the homogeneous solution of Eqs. ([3]). How- 
ever, as shown in Fig. [1] the entropy per site in the su- 
percooled liquid becomes negative when the density is 
increased above ps=o — 0.1757 (or the chemical poten- 
tial is increased above ps=o — 7.07). As a consequence, a 
thermodynamic phase transition must occur at a density 
p < ps=o- In fact, in the mean-field approximation we 
find that the system undergoes a phase transition toward 
a 1RSB glassy phase, which can be analyzed by taking 
into account the existence of many different local minima 
(or pure states) of the free energy. Since in this case the 
local fields can fluctuate from pure state to pure state, 
this situation is described by a probability distribution 
P(h, v) that the fields hi and Vi on the site i equal h 
and v for a randomly chosen pure state. Using the cav- 
ity method we find that this function satisfies the 
self-consistent equation: 



From these fields one can obtain the free energy, j3F = 
— \nZ, as a sum of site and link contributions [14j : F = 
AF S - (k + l)AFi/2. The contribution from the bond 
between two branches with root sites i\ and 12 is: 



e -0AF l _ 1 + e p(h il +ih 2 ) + e 0(h i2 +v n )^ 



while the contribution from the addition of a site i con- 
nected with k + 1 branches with root sites ij reads: 



fe+i 



-0AF s 



(5) 



3 = 1 



3=1 

(6) 

where Af is a normalization constant, hi and m are the lo- 
cal cavity fields obtained when merging k branches which 
carry the fields (hi^v^) [via Eq. ([3])], and m 6 [0,1] is 
a Lagrange multiplier which turns out to be the usual 
1RSB parameter, fixed b y th e maximization of the free 
energy with respect to it [Til [H| . 

The 1RSB cavity equation, Eq. ©, can be solved ex- 
actly in the close packing limit (/3 — > oo), where the re- 
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FIG. 2: Schematic drawing of the location of the rele- 
vant densities emerging from the mean-field solution: melt- 
ing density, p m ~ 0.1364, dynamical transition density, pd ~ 
0.1688, Kauzmann transition density, pK — 0.1750, den- 
sity at which the entropy of the supercooled liquid vanishes, 
ps=o — 0.1757(~ Pk), maximum density of the crystalline 
packing, p m ax = 0.2. The figure also show that the configu- 
rational entropy, E, defined as the logarithm of the number 
of pure states, is positive for p d < p < pK- 



cursion relations, Eq. ([3]), simplify to: 

hj = 1 +} j v ij - V(h i:j ), v 



-V{h H ) (7) 



with V(hi ) = max{/ii }0(max{/ii }) [we have set fi = 
1]. These equations yield an exact ansatz for the cavity 
fields probability distribution: P{h,v) = Y^1=iPr3(h — 

h r )5(v - V r ), With /ll,2,3,4 = 1,0,-1,-2, Ul,2,3,4 = 0, 

/i5,6,7 = 0,-1, -2, and w 5j6:7 = -1. By defining P = pi, 
R = P2 + P3 + Pi-, and Q = p 5 + p e + Pi = 1 - P - R, 
Eq. ([6]) becomes: 



P 



R 3 [e y - (e v - 1)(1 - P) 3 ]" 



(8) 



R = [(1 - P) 3 - i? 3 ] [e* - {e y - 1)(1 - P) 3 \ 



where y = lim/^oo (3m. In terms of P and P, the 1RSB 
link and site contribution to the free energy read: 

A(/)i[y} = -y- 1 ln{l + (e y - l)P 2 + 2e v PR} (9) 
AcPM = -y- 1 ln{e^ + (e^-l)[i? 4 -(l-P) 4 ]}. 

The free energy is then given by <f>[y] — A<fi s [y] — 
2<f>i[y], from which we can compute the complexity £ = 
y 2 d(j)[y]/dy and the density p = d{ycj)[y\) / dy . 

The finite pi solution of the 1RSB cavity equations can 
be found numerically, using the population dynamics al- 
gorithm 14 1 . In agreement with the results of Refs. H,[l], 
the mean-field solution of the model exhibits an ideal 
glass transition of the same nature of that found in mean- 
field discontinuous spin glasses [ll[ . We first find a purely 
dynamical transition at pd — 0.1688 {fid — 5.86), where 
a non trivial solution of Eq. ^ appears for the first 
time, signaling the emergence of an extensive number of 
metastable states (which, in mean-field, trap the dynam- 
ics for an infinite time). A solution of the 1RSB equation 
becomes thermodynamically stable at a higher density, 
px — 0.1750 {fiK — 6.92), where a thermodynamic tran- 
sition to a 1RSB glassy phase takes place. The relevant 
densities emerging from the mean-field approximation are 
reported in Fig. [21 showing that pk is strikingly close to 
Ps=o- 
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FIG. 3: Expansion around the mean-field solution. Starting 
from the mean-field limit {L = 1), bigger finite-dimensional 
cells of size L = 2 and L = 3 are considered. 



X 



X 



&<xr 



X 
X 



X 
XX 



FIG. 4: Case L = 2: the figure depicts a configuration in 
which a particle (red circle) occupies one of the two forward 
sites of the branch ending on the cell i, i.e., one of the two 
sites which are not on the edge of the cell where the link is 
missing. The sites whose distance from the particle is less 
or equal than two lattice spacings cannot be occupied (these 
sites are marked by crosses). 



In the following we consider a systematic expansion 
around the mean-field theory, which takes into account 
the actual structure of the two dimensional square lattice. 
More precisely, we consider cells of size L = 2, 3 within 
which the model is solved "exactly" and we use those 
cells as vertex of the mean-field theory on the random 
graph, as depicted in Fig. [3] (see also Refs. [H, [III). In 
the limit L — > oo the exact solution of the model should 
be achieved. This method allows to include in an exact 
fashion short range spatial correlations, as correction to 
the mean- field limit. Since in glassy systems there is 
no diverging equilibrium length scale, this expansion is 
expected to be reliable and effective. 

In order to solve the L — 2 case, let us consider a 
branch of the tree ending on the 2x2 square cell i, and 
denote by ij, j € {1, . . . , k}, the roots of the branches 
connected to the cell i, as shown in Fig. 2J We call 
(resp., flW) the partition function of the branch 
restricted to the configurations in which the cell is oc- 
cupied by a particle in one of the two "forward" (resp., 
"backward") sites, i.e., one of the two sites which are not 
on the edge of the cell where the link is missing (resp., 
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one of the two sites which are on the edge of the cell 

where the link is missing). We also define U as the 
partition function of the branch restricted to the config- 
urations in which the cell is occupied by a particle in 
one of the two forward sites, with its "left" (or, equiva- 
lently, "right") neighbor cell constrained to be not occu- 
pied in one of the two backward sites. Finally, we need 
to introduce , defined as the partition function of the 
branch restricted to configurations in which the cell is 

empty, and Zq , corresponding to the partition func- 
tion of the branch restricted to configurations in which 
the cell is empty with its "left" (or, equivalently, "right") 
neighbor cell constrained to be not occupied in one of the 
two backward sites. Introducing four local cavity fields, 

defined as (3 Ui = ]n(U® /zfi), /3dj = \n{D^/Z^), 

0Vi = ln(I7 (i) /Zft), and (3 Pi = \n(Z~ Q {i) JZ^), within the 
RS homogeneous ansatz for the supercooled liquid, the 
following algebraic recursion relations are derived: 

e U = eP^WZ- 1 , = e^T 2 YZ~ 1 

e t3d = e J3n T ( e Pd + W ^YZ- 1 (10) 

e 0p = [ e W(We 0d + Y ■F)+Y{e pd + F)W]Z-\ 
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FIG. 5: Entropy per site, S, as a function of the density, p. 
L increases along the direction of the arrow: L — 1 (continu- 
ous curve), L = 2 (dashed curve) and L = 3 (dotted-dashed 
curve). The density at which the entropy of the supercooled 
liquid vanishes, ps=o, increases as L is increased. 



where T 



+ e 



0V 



Y = 1 + 2e /3u , W = 
1 + eJ M + 2e /3u , F = W + e 0d , and Z = 
[e 0d (We 1311 + Y ■ F) + W 2 (e^ + F)] . In terms of the 
local cavity fields, the contribution to the free energy 
from the addition of a link between two cells, AFi, and 
from the addition of a cell, AF S , can be computed: 



-0AFi _ 



e ' = Y z + 4e 0d (e Pv + er / 



(11) 

2e l3(d+v \ 



from which one can derive the free energy per site F = 
(AF S — 2AF))/4, the density of particles p = (rij) = 
— (<9F/9/z)/4, and the entropy per site S — —(3F — (3pp. 

The recursion relations for the case L — 3 can be de- 
termined using a similar procedure. In Fig. [5J the en- 
tropy, S, is plotted as a function of the density, p, for 
L = 1, L = 2 and L = 3, in the homogeneous solution 
for the supercooled liquid. The figure shows that the ex- 
pansion around the mean-field limit systematically mod- 
ifies the results. In particular, we note that the density 
at which the entropy of the supercooled liquid vanishes, 
Ps=o(~ Pk), moves toward higher densities when L is 
increased. As a consequence, the instability of the super- 
cooled liquid (and, therefore, the transition to the glassy 
phase) is displaced toward the maximum density when 
bigger two dimensional cells are considered. In fact, we 
find that ps=o approaches very nicely p m ax = 0.2 as a 
power law: p s =o{ L ) - O- 2 - 0.024 L~ - 44 . 

Further insights can be gained by studying the value of 
the entropy of the supercooled liquid in the zero tempera- 



ture limit, defined as S Q 



lira,. 



S . S 



is negative in 

the mean-field approximation [Soo(L = 1) ~ —0.291] and 
increases systematically as L is increased. We find that 



Sao nicely approaches zero as a power law as a function 
of L: ~ —0.291 L~ 36 . 

These results clearly hints that the instability of the 
supercooled liquid and the thermodynamic transition to 
the glassy phase is smeared out in two dimensions, as 
one includes corrections to the mean-field theory: the 
entropy of the supercooled liquid seems to vanish only 
at the maximum density of the crystalline state. This 
is in agreement with the recent findings of Refs. [H, 6|. 
Our results are also in agreement with those of Ref . 7| , 
where, by employing a suitable Monte Carlo algorithm, 
the authors show that there is no evidence for a ther- 
modynamic phase transition up to very high densities in 
two dimensional (polidispcrse) hard disks; the glass is 
thus indistinguishable from the liquid on purely thermo- 
dynamic grounds. Note, however, that in Ref. [T(| the 
authors state that numerical claims in favour of and/or 
against a thermodynamic glass transition must be consid- 
ered carefully, due to the difficulties to probe the system 
close enough to px- 

If the liquid has to be a good solution in the L — > 
oo limit up to p m ax, its pressure must diverge at this 
point. In fact, we find that for every value of L, the 
pressure of the supercooled liquid diverges for p — > oo, 
i.e., the entropy of the liquid approaches Soo at p max 
with a vertical slope. Since, as discussed above, is 
extrapolated to zero and ps=o is extrapolated to p ma x as 
L is increased, it seems reasonable to expect that in the 
large L limit the entropy of the liquid vanishes at p max 
and that, consistently, the pressure diverges at this point. 

However, it is important to highlight that this is not a 
proof but just a hint of the absence of the glass transition 
in hard disks: as a matter of fact, ps=o is only an upper 



5 



bound to pk- Thus it might be possible that, in the limit 
of large L, p K < p max even if p s=0 -> p max . In princi- 
ple, one should compute pk for different values of L, by 
solving the 1RSB equations, which is, unfortunately, a 
hard numerical task. Nevertheless, given the closeness of 
Pk to ps=o for L = 1, and the consistency of the liquid 
solution when extrapolated for large L up to p ma x, one 
might guess that the possibility described above is un- 
likely and that ps=a provides a good estimation for px 
also for bigger L. 

In conclusion, we have presented an analytical study 
of a system of identical hard spheres, focusing on the 
case of hard disks on a square lattice. The mean-field 
version on the model exhibits an ideal glass transition 
of the same kind of that found in mean-field discontinu- 
ous spin glasses Q. Nevertheless, by considering a sys- 
tematic expansion around the mean-field solution able 
to take into account short range correlations in an exact 
fashion, we have shown that such glass transition seems 
to be smeared out in two dimensions, confirming the re- 
sults of Refs. [E S 0] ■ Note that the results presented 
here are also in agreement with the recently discovered 
mapping of glass forming system to Ising spin glasses 
in an external magnetic field [13], according to which 



there should not be a thermodynamic glass transition in 
dimensions less than six. Since there cannot be a dy- 
namical glass transition without a thermodynamical one 
(provided that the dynamics satisfies the detailed bal- 
ance), we finally argue that there is no structural arrest 
in two dimensional hard disks at a density smaller than 
the one of the crystalline packing. This analysis could 
be generalized to the three dimensional case (which is 
the most relevant for supercooled liquids) and to higher 
spatial dimensions, providing a direct way to investigate 
the lower critical dimensionality of the mean-field ideal 
glass scenario. This study could prune down the number 
of candidate theories for the glass transition. 
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